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Abstract 

The  present  report  investigates  using  simultaneous  observations  of  coincident  gravity  waves  and  sprites  to  establish  an  upper 
limit  on  sprite-associated  thermal  energy  deposition  in  the  mesosphere.  The  University  of  Alaska  operated  a  variety  of  optical 
imagers  and  photometers  at  two  ground  sites  in  support  of  the  NASA  Sprites99  balloon  campaign.  One  site  was  atop  a  US 
Forest  Service  lookout  tower  on  Bear  Mt.  in  the  Black  Hills,  in  western  South  Dakota.  On  the  night  of  18  August  1999  we 
obtained  from  this  site  simultaneous  images  of  sprites  and  OH  airglow  modulated  by  gravity  waves  emanating  from  a  very 
active  sprite  producing  thunderstorm  over  Nebraska,  to  the  Southeast  of  Bear  Mt.  Using  25  s  exposures  with  a  bare  CCD 
camera  equipped  with  a  red  filter,  we  were  able  to  coincidentally  record  both  short  duration  ( <  10  ms)  but  bright  ( >3  MR) 
N2  IPG  red  emissions  from  sprites  and  much  weaker  (~1  kR),  but  persistent,  OH  Meinel  nightglow  emissions.  A  time 
lapse  movie  created  from  images  revealed  short  period,  complete  360°  concentric  wave  structures  emanating  radially  outward 
from  a  central  excitation  region  directly  above  the  storm.  During  the  initial  stages  of  the  storm  outwardly  expanding  waves 
possessed  a  period  of  t«10  min  and  wavelength  2%50  km.  Over  a  1  h  interval  the  waves  gradually  changed  to  longer  period 
t«11  min  and  shorter  wavelength  a^40  km.  Over  the  full  2  h  observation  time,  about  two  dozen  bright  sprites  generated 
by  the  underlying  thunderstorm  were  recorded  near  the  center  of  the  outwardly  radiating  gravity  wave  pattern.  No  distinctive 
OH  brightness  signatures  uniquely  associated  with  the  sprites  were  detected  at  the  level  of  2%  of  the  ambient  background 
brightness,  establishing  an  associated  upper  limit  of  approximately  AT  <  0.5  K  for  a  neutral  temperature  perturbation  over 
the  volume  of  the  sprites.  The  corresponding  total  thermal  energy  deposited  by  the  sprite  is  bounded  by  these  measurements 
to  be  less  than  ~1  GJ.  This  value  is  well  above  the  total  energy  deposited  into  the  medium  by  the  sprite,  estimated  by  several 
independent  methods  to  be  on  the  order  of  ~1-10  MJ. 

©  2003  Elsevier  Science  Ltd.  All  rights  reserved. 
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1.  Introduction 

Sprites  are  transient  (several  to  10’s  of  ms)  optical 
emissions  generated  within  the  mesosphere  (50-90  km) 
by  large  lightning  discharges  (e.g.,  Sentman  and  Wescott, 
1993;  Lyons,  1994;  Fukunishi  et  al.,  1996;  cf.  Sentman, 


"Corresponding  author.  Tel.:  +1-907-474-6442;  fax:  +  1-907- 
474-7290. 

E-mail  address :  dsentman@gi.alaska.edu  (D.D.  Sentman). 


1998).  Several  mechanisms  have  been  proposed  for  sprite 
production  (see  review  by  Rowland,  1998),  all  based  on  im¬ 
pact  excitation  of  ambient  molecular  species  by  electrons 
accelerated  in  transient  mesospheric  electric  fields  following 
a  lightning  discharge  in  the  underlying  thunderstorm.  The 
most  generally  accepted  mechanism  is  the  quasi-electrostatic 
mechanism  of  Pasko  et  al.  ( 1997a).  Imaging  and  photomet¬ 
ric  measurements  reveal  a  continuous  range  of  sprite  bright¬ 
nesses,  from  a  few  tens  of  kR  for  small  events  to  several 
MR  for  the  brightest  forms  (Moudry  et  al.,  2003).  Studies  of 
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sprite  occurrence  over  thunderstorms  in  the  central  US  show 
that  sprites  tend  to  occur  predominantly  over  the  trailing 
stratoform  regions  of  large  mesoscale  convective  systems 
(MCSs)  of  size  greater  than  ~20,000  km2  (Lyons,  1996), 
although  studies  in  other  regions  of  the  world  show  them 
also  to  occur  over  smaller  storms  (Fukunishi  et  al.,  1999; 
Su  et  al.,  2002). 

Spectroscopic  studies  have  shown  the  optical  emissions 
of  sprites  are  predominantly  from  the  red  ^B^rig^v')  — > 
A3Zu(v")  IPG  (First  Positive)  bands  of  neutral  molecular 
nitrogen  (Mende  et  al.,  1995;  Hampton  et  al.,  1996;  Heavner 
et  al.,  2000).  Weak  ion  signatures  have  been  reported  in  both 
red  Nj(/42/7  — ►  X2Zg)  Meinel  emissions  (Bucsela  et  al., 
1998;  Morrill  et  al.,  1998)  and  blue  N2+(£2r+  ->  X2Zg  ) 
1NG  (First  Negative)  emissions  (Armstrong  et  al.,  1998; 
Suszcynsky  et  al.,  1998;  Takahashi  et  al.,  1998;  Armstrong 
et  al.,  2000).  Electron  impact  excitation  and  ionization  of 
O2  are  also  believed  to  occur  (Pasko  et  al.,  1997a),  but  col- 
lisional  quenching  suppresses  optical  emissions  from  these 
species.  The  characteristic  electron  energy  of  a  few  eV  lead¬ 
ing  to  the  N2  IPG  emissions  also  leads  to  excitation  of  ro¬ 
tational  and  vibrational  states  of  the  ambient  species,  within 
both  the  observed  transition  states  of  N2  (Cartwright,  1978; 
Green  et  al.,  1996)  and  a  larger  group  of  quenched  states  of 
both  N2  and  O2  (Morrill  et  al.,  2002;  Bucsela  et  al.,  2003). 
Thus,  while  the  N2  IPG  emissions  constitute  the  bulk  of  the 
optical  energy  observed,  they  are  likely  to  represent  only 
a  small  fraction  of  the  total  energy  deposited  in  the  meso¬ 
sphere  at  the  location  of  a  sprite. 

Establishing  the  total  energy  deposited  by  sprites  in  the 
middle  and  upper  atmosphere  is  important  within  the  context 
of  the  overall  energy  budget  of  this  region,  including  inputs 
from  other  external  sources,  such  as  by  solar  EUV  or  au¬ 
roral  particle  precipitation  from  above,  or  by  gravity  waves 
from  below.  Stenbaek-Nielsen  et  al.  (2000)  have  presented 
observations  suggesting  local  composition  effects  in  sprite 
tendril  formation,  but  the  relative  importance  of  large  scale 
composition  or  dynamical  effects  induced  in  the  embedding 
medium  depends  significantly  on  the  total  amount  of  energy 
available  compared  to  other  sources.  In  situ  diagnosis  of  the 
energy  deposition  associated  with  the  electrical  perturbation 
of  sprites  is  hampered  by  the  difficulty  of  placing  instru¬ 
ments  within  the  discharge  region  at  the  time  of  the  events. 
One  is  thus  led  to  consider  signatures  of  secondary  effects 
produced  in  the  medium  by  the  electrical  perturbation  that 
produces  sprites,  such  as  the  ambient  airglow  layers  in  the 
upper  atmosphere. 

Sprites  typically  extend  to  altitudes  of  85-90  km  (Sent¬ 
man  et  al.,  1995),  and  thus  penetrate  into  the  airglow  re¬ 
gions  maintained  in  the  mesosphere  by  a  balance  of  vari¬ 
ous  chemical  and  transport  processes.  The  visible  volume  of 
large  sprites  often  exceeds  104  km3  for  large  events  (Sent¬ 
man  and  Wescott,  1993).  Their  duration  is  short  (few  ms), 
so  thermal  energy  deposited  in  the  neutral  atmosphere  within 
a  sprite  would  be  expected  to  produce  an  impulsive  pres¬ 
sure  pulse  that  propagates  laterally  outward  as  an  acous¬ 


tic  wave  or  gravity  wave.  Pasko  et  al.  (1998)  have  esti¬ 
mated  the  neutral  heating  perturbation  within  streamer  struc¬ 
tures  to  be  AT/T~0.2-2%  at  altitudes  at  altitudes  of  50 
-60  km.  High  resolution  imaging  of  sprite  structures  have 
revealed  them  to  be  composed  of  multiple  thin  stream¬ 
ers  (Gerken  et  al.,  2000),  so  the  “filling  factor”,  or  the 
fractional  volume  occupied  by  streamers  within  the  body 
of  a  sprite,  is  very  small,  estimated  to  be  ~10“4-10-3. 
The  associated  neutral  temperature  perturbation  averaged 
over  the  full  body  of  the  sprite  is  then  (0.3-3  K)  x  (10“4- 
10“ 3 )~3  x  10~5  —  3  x  10-2  K.  For  a  sprite  volume  of 
1012  m3  and  neutral  density  of  2  x  1020/m3  at  85  km,  this 
corresponds  to  a  total  energy  range  ~0. 1-1 0  MJ.  Other  work 
based  on  scaling  of  observed  optical  emission  line  ratios 
has  suggested  that  the  total  energy  deposited  in  the  meso¬ 
sphere  by  sprites  could  be  as  large  as  ~1  GJ  (Heavner 
et  al.,  2000). 

Given  the  large  range  of  uncertainty  in  the  energy  depo¬ 
sition,  spanning  some  4  orders  of  magnitude,  it  was  deemed 
worthwhile  to  search  for  a  possible  temperature  perturba¬ 
tion  signature  by  way  of  correlated  effects  on  the  nightglow 
emissions  (Sentman  et  al.,  1999),  both  as  a  way  to  narrow 
the  range  of  uncertainty  and  to  clarify  the  nature  of  possible 
secondary  effects  on  airglow  processes.  Positive  detection  of 
sprite-correlated  nightglow  pulses  would  provide  a  method 
for  indirectly  determining  the  underlying  thermal  energy  de¬ 
position  associated  with  the  sprites,  as  well  as  other  related 
processes,  while  a  negative  result  would  permit  establishing 
an  upper  limit  on  the  energy  deposition. 

Rowland  et  al.  (1996)  and  Pasko  et  al.  (1997a,  b)  have 
considered  the  converse  problem,  that  of  the  manner  in 
which  effects  of  gravity  waves  on  the  background  might 
play  a  role  in  determining  the  ignition  location  or  structure 
of  sprites.  Siefring  et  al.  ( 1 998 )  have  presented  examples  of 
apparent  close  association  between  sprites  and  short  period 
gravity  wave  structures  observed  during  one  of  the  EXL98 
aircraft  flights,  suggesting  a  possible  linkage  between  the 
two. 

In  Section  2  the  linkage  between  gravity  waves  and  OH 
emissions  is  reviewed.  Section  3  describes  the  experiment 
that  yielded  the  observations,  which  are  presented  in  Section 
4.  Section  5  describes  the  analysis  of  the  observations,  their 
interpretation  in  terms  OH  Meinel  emissions,  and  the  cor¬ 
responding  energy  limits  on  sprite  deposition  derived  from 
these  observations.  A  short  discussion  on  related  issues  is 
given  in  Section  6.  Results  and  some  additional  comments 
about  the  observations  conclude  the  paper  in  Section  7. 

2.  OH  chemistry  and  wave  dynamics  background 

Of  particular  interest  to  this  study  is  the  OH  nightglow 
layer  located  at  ~85  km,  first  reported  by  Meinel  (1950). 
The  production  of  vibrationally  excited  hydroxyl  OH(t  )  in 
the  nighttime  mesosphere  is  due  to  a  set  of  coupled  chemi¬ 
luminescent  reactions  whose  basic  processes  involve  ozone 
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destruction  by  atomic  hydrogen  (Bates  and  Nicolet,  1950). 
The  set  of  reactions  needed  to  explain  observations  has  sub¬ 
sequently  been  expanded  to  include  other  odd  oxygen  and 
odd  hydrogen  species  (Makhlouf  et  al.,  1995), 

H  +  O3  — ►  OH(t;  ^  9)  +  02, 

*,  =  1.4  x  10-,0e-470/r. 

0  +  H02  —  OH(t>  s;  6)  4-  02, 

*2  =  3  x  10_", 

O  +  O2  4-  M  —>03+  Af, 

*3  =  6  x  10'34(300/r)2J, 

H  +  02  +  A/  —  H02+M, 

*4  =  5.7  x  10"32(300/r)16, 

OH(u)  +  O  — ►  H  +  02 

k5  =  03(f)  x  10-11,  as(v  =  0)  =  3.9, 

as(v  =  1)  =  10.5,  as( v  >  1 )  =  25,  (1) 

where  T  is  the  ambient  neutral  temperature  in  K ,  and  k\ ...  £5 
are  rate  coefficients  in  cm3/s  for  two-body  reactions,  and 
cm6/s  for  three-body  reactions  involving  a  third  species  A/, 
either  N2  or  O2.  For  the  first  two  reactions  the  total  rate 
coefficient  for  production  of  all  vibrational  states  v  is  listed. 
Near  the  mesopause  the  neutral  temperature  is  T~170  K. 
The  vibrational ly  excited  states  can  be  quenched  by  collision 
with  major  species  or  chemically  by  the  OH(t;)  -1-  O  reaction 
above,  but  the  major  loss  over  most  of  the  region  of  interest 
is  by  radiation, 

OH(z;)  OH(t-  -  At;)  +  hi\  Av>r- ap,  (1') 

where  Ar,c-&v  is  the  Einstein  coefficient  for  spontaneous 
emission,  v  =  1,...,9,  and  0<Av^v.  The  observed 
radiance  from  the  vibrational  level  v  to  level 

v-Av  transition  is  then  the  integral  of  the  volume  emission 
rate  i4*,p-Ar[OH(t;)]  along  the  line-of-sight,  where  [OH(t;)] 
is  the  density  of  species  OH  in  vibrational  state  v.  The 
nighttime  airglow,  or  nightglow,  generated  by  these  reac¬ 
tions  produces  bright  emissions  (hundreds  of  kR)  from  the 
OH(X2n,  At;  =  2)  first-overtone  Meinel  sequences  at  wave¬ 
lengths  1.4 -2.2  pm,  with  the  total  intensity  of  all  emission 
sequences  exceeding  several  MR  (Chamberlain,  1961). 

Gravity  waves  produced  by  convective  core  motions  of 
vigorous  thunderstorms  propagate  into  the  high  atmosphere 
and  constitute  a  major  source  of  mechanical  energy  and 
momentum  into  this  region  from  below  (Swenson  and  Liu, 
1998).  Among  other  effects,  these  thunderstorm-generated 
gravity  waves  can  modulate  the  optical  emissions  in  the 
nightglow  layer  (Krassovsky,  1972)  that  are  detectable  from 
the  ground  (e.g.  Taylor  and  Hapgood,  1988;  Taylor  and  Hill, 
1991;  Turnbull  and  Lowe,  1991 )  and  from  orbiting  imaging 
platforms  (Dewan  et  al.,  1998;  Picard  et  al.,  1998).  The 


emission  intensity  of  the  OH  nightglow  is  sensitively 
dependent  on  small  wave-induced  temperature  fluctua¬ 
tions  (Makhlouf  et  al.,  1995)  and  readily  responds  to 
thunderstorm-generated  gravity  wave  temperature  fluctua¬ 
tions  to  produce  brightness  variations  such  as  reported  by 
Taylor  et  al.  (1995). 

The  effect  of  a  gravity  wave  on  the  OH  nightglow  emis¬ 
sions  is  determined  by  computing  the  linear  response  of 
the  coupled  set  of  Eqs.  ( 1 )  and  ( l')  to  a  small  temperature 
perturbation  AT.  Because  the  rate  coefficients  in  the  reac¬ 


tions  above  possess  different  temperature  dependencies,  the 
steady-state  emissivity  is  also  temperature  dependent.  The 
differential  response  of  the  OH  Meinel  Emission  intensities 
low  to  small  temperature  perturbations  was  first  studied  by 
Krassovsky  ( 1 972 )  and  may  be  expressed  in  terms  of  the 
differential  temperature  perturbation  as 


d/OH  dT 
/oh  ***  T  ' 


(2) 


where  rj  is  a  factor  that  depends  on  wave  parameters  in¬ 
cluding  frequency  and  phase  speed.  For  mesosphere  condi¬ 
tions  and  the  set  of  reactions  in  Eqs.  (1)  and  (1')  above, 
the  factor  rj  ranges  from  \rj\  =  10  for  cph  =  40  m/s  to  \rj\  = 
5  for  cph  =  140  m/s  (Makhlouf  et  al.,  1995).  Hence,  a 
small  temperature  change  produces  a  correlated  OH  optical 
emissivity  response  amplified  by  a  factor  of  5-10.  To  scale 
optical  emission  variations  to  the  underlying  temperature  ef¬ 
fects  requires  knowledge  of  the  amplitude  and  phase  of  the 
rj  parameter,  which  can  be  determined  from  the  measurable 
parameters  of  the  wave  modulating  the  OH  emission. 

In  the  present  study  we  searched  for  transient,  in  situ  heat¬ 
ing  effects  at  the  mesopause  produced  by  sprites  through 
their  effects  on  the  modulation  of  the  nocturnal  OH  airglow 
emissions,  similar  to  the  effects  produced  by  gravity  wave 
modulation  of  the  OH  emissions.  While  making  the  mea¬ 
surements  of  sprites  above  a  very  vigorous  thunderstorm  we 
unexpectedly  captured  a  well-defined  example  of  an  out¬ 
wardly  expanding,  gravity  wave  driven  concentric  ripple 
pattern  in  OH  emissions  driven  by  the  thunderstorm,  simi¬ 
lar  to  effects  first  reported  by  Taylor  and  Hapgood  (1988), 
and  further  described  by  Dewan  et  al.  (1998).  The  wave¬ 
length  and  period  extracted  from  these  narrow  band  waves 
provided  a  key  parameter  for  characterizing  the  perturba¬ 
tion  response  of  OH  needed  to  evaluate  the  sprite  heating 
parameters. 


3.  Experiment 

Observations  were  performed  during  August  1999  as 
part  of  the  NASA  Sprites99  campaign  (Bering  et  al.,  1999; 
Bering  et  al.,  2002)  at  a  US  Forest  Service  fire  observation 
tower  on  Bear  Mt.  (latitude  43.88N,  longitude  103.75W, 
altitude  2150  m),  located  in  the  Black  Hills  near  the  city 
of  Custer,  South  Dakota.  The  observing  camera  was  a 
Photometries  SenSys  Model  400  768  x  512  pixel  x  12  bit 
unintensified  CCD  system  equipped  with  a  Kinoptic 
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CCD  Camera  and  Red  Filter  Responses 


Fig.  1 .  CCD  and  cutoff  filter  responses  as  functions  of  wavelength, 
and  the  combined  response  pass  band.  Also  shown  are  locations 
of  the  band  origins  of  the  OH  (X2/7)  Av  =  4,5  Meinel  sequences 
within  the  combined  response  pass  band. 


5.7  mm  //1.8  TEGEA  wide  angle  lens.  The  CCD-lens 
combination  provided  a  60° H  x  40° V  degree  field  of 
view.  We  used  a  Schott  RG-715  715-nm  red  filter  to  pass 
a  selected  band  of  wavelengths  in  the  NIR  region  of  the 
spectrum.  When  combined  with  the  camera  response,  the 
camera-filter  combination  (NIR  CCD)  possessed  an  overall 
optical  pass  band  of  715-920  nm  FWHM.  Fig.  1  shows  the 
respective  response  functions  of  the  camera-filter  system. 

The  OH  nightgiow  emissions  we  observed  are  not  the 
brightest  among  the  several  types  of  nocturnal  skyglow. 
As  noted  above,  the  brightest  emissions  occur  in  the 
short-wave  infrared  (SWIR)  (X2n,Av  =  2)  first-overtone 
Meinel  sequences  1. 4-2.2  pm.  However,  the  OH  intensities 
are  significant  (several  tens  kR)  in  the  shorter  wavelength 
700-900  nm  near-infrared  (NIR)  OH(X2n,Av  =  3,4,5) 
Meinel  sequences  (Yee  et  al.,  1991 )  and  are  the  emissions 
reported  here.  The  CCD  cooling  requirements  for  detection 
of  the  NIR  emissions  are  much  less  stringent  than  for  SWIR, 
and  make  it  possible  to  use  standard  CCD  technology 
operating  near  the  upper  limits  of  the  camera  wavelength 
sensitivity.  Thus,  one  may  use  bare  (unintensified)  CCD 
technology  for  sufficiently  long  exposure  times  (tens  of 
seconds)  (e.g.,  Taylor  and  Hill,  1991;  Taylor  et  al.,  1991a, 
b)  with  a  relatively  modest  (T~0°C)  amount  of  cooling 
and  is  the  method  we  used  for  the  present  observations. 

Table  1  lists  the  principal  OH  Meinel  bands  and  their  loca¬ 
tions  within  the  overall  instrument  pass  band  700-900  nm. 
The  pass  band  included  the  N2(B3ng  — >  A3 If) Av  =  2 
(740—780  nm)  and  Av  =  1  (850-900  nm)  sequences,  so 


Table  1 

OH  ( X2n )  Meinel  vibration-rotation  bands  720 <2 <  1000  nm 


Band 

k  (nm) 

(Q-branch) 

Relative 

intensity 

At;  =  3 
(3-0) 

966 

2 

Av  =  4 
(9-5) 

(R-branch) 

1000 

3 

(8-4) 

940 

8 

(7-3) 

883 

7 

(6-2) 

834 

5 

(5-1) 

792 

6 

(4-0) 

753 

3 

At;  =  5 
(9-4) 

776 

10 

(8-3) 

728 

6 

Source :  Pearse  and  Gaydon  (1976). 


the  system  was  also  sensitive  to  the  sprite  IPG  emissions. 
Hence,  we  were  able  to  image  both  OH  nightgiow  emissions 
and  sprites  simultaneously. 

The  operational  mode  of  the  present  experiment  involved 
performing  25  s  exposures  followed  by  a  5  s  period  to  trans¬ 
fer  the  digital  images  from  the  camera  to  hard  disk,  yielding 
an  effective  image  frame  rate  of  1/30  fps.  The  camera  was 
connected  via  an  interface  card  to  a  desktop  computer.  Im¬ 
age  acquisition  and  readout  were  controlled  with  a  custom 
script  running  within  the  main  Photometries  camera  control 
program.  The  computer  system  timer,  updated  continuously 
over  the  Internet  to  maintain  100  ms  accuracy  absolute,  was 
used  to  synchronize  image  acquisition.  The  resultant  set  of 
images  was  processed  off  line  for  histogram  equalization  us¬ 
ing  an  average  histogram  for  the  entire  sequence  as  a  refer¬ 
ence,  and  combined  into  a  digital  video  clip  for  further  anal¬ 
ysis  and  presentation.  A  separate  intensified  CCD  (TCCD) 
television  system  was  also  operated  during  the  campaign  to 
obtain  images  of  many  of  the  brighter  sprites  that  registered 
in  the  NIR  images. 

4.  Observations 

On  the  night  of  18  August  1999  over  a  2  h  interval  we 
obtained  very  clear  images  of  both  gravity  wave  modu¬ 
lated  nightgiow  emissions  and  sprites  associated  with  a  very 
intense  thunderstorm  over  Nebraska.  The  thunderstorm  is 
shown  in  GOES  infrared  imagery  with  overlaid  lightning 
strike  information  in  Fig.  2.  Over  the  3  h  period  0400-0700 
UT,  the  center  of  storm  convection  and  lightning  activity 
moved  approximately  200  km  from  south-central  Nebraska 
to  eastern  Nebraska/Iowa  at  an  average  speed  of  about  60 
-70  km/h  (17-19  m/s).  This  storm  track  was  roughly  or¬ 
thogonal  to  the  observing  line  of  sight,  which  approximately 
preserved  the  distance  to  the  storm  during  the  observing 
interval.  Observing  conditions  were  close  to  ideal,  with  a 
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Fig.  2.  GOES-8  10-12  pm  temperature  maps  of  a  thunderstorm  over  central  Nebraska-western  Iowa  and  associated  lightning  studied  during 
the  night  of  18  August  1999.  The  two  maps  show  the  same  storm  separated  by  three  hours.  Superimposed  atop  the  maps  are  5  min  averages 
of  lightning  activity,  where  blue  denotes  negative  CG  (cloud-to-ground)  lightning,  and  red  denotes  positive  CG  lightning.  Adapted  from 
near-real  time  maps  downloaded  from  http://www.weather.com. 


Clouc 
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Looking  ESE  from  Bear  Mt  SD 
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Photometries  SenSys  400 
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720  nm  red  filter 


Cloud  Bank 


Fig.  3.  Outwardly  expanding  concentric  ripples  observed  in  OH  emissions  using  25  s  time  exposure  in  the  NIR.  The  field  of  view  is 
approximately  60  H  x  40  V  deg.  The  dashed  lines  labeled  “analysis  strip”  delineate  the  region  used  to  create  Fig.  5.  The  city  lights  of 
Custer,  SD  and  lights  from  the  Crazy  Horse  Monument  are  labeled  at  the  bottom  of  the  large  plot.  The  stars  in  the  figure  are  used  to  fix 
the  viewing  direction  using  the  Smithsonian  Star  Catalog.  The  sprite  on  the  main  image  is  shown  in  the  inset  in  more  detail  (recorded  by  a 
separate  ICCD  camera). 


clear  sky  above  the  observing  site  permitting  visual  access 
to  nightglow  and  sprite  activity  in  the  mesosphere  above  the 
thunderstorm. 

Fig.  3  shows  one  of  the  NIR  CCD  images  obtained  during 
the  observation  interval.  In  the  foreground  near  the  bottom  of 
the  large  image  are  seen  the  lights  of  the  city  of  Custer,  SD, 


and  to  the  left  the  display  lights  of  the  Crazy  Horse  Monu¬ 
ment.  Left  lower  center  there  is  a  large  sprite,  shown  in  high 
resolution  as  a  cutout  from  the  ICCD  TV  camera.  The  stars 
are  clearly  visible  at  the  top  of  the  figure.  The  hazy  ellipti¬ 
cal  features  in  the  lower  central  portion  of  the  figure  labeled 
as  “Concentric  Expanding  Gravity  Wave  Ripples”  are  the 
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Fig.  4.  GOES-8  weather  image  with  image  of  Fig.  3  mapped  to 
the  85  km  tangent  surface  and  overlaid  on  top  of  it.  The  outward 
expanding  ripples  are  seen  to  be  nearly  centered  on  the  underlying 
thunderstorm.  See  text  for  description  of  the  mapping  projection 
used. 


gravity  wave  signatures  we  seek.  When  the  30  s  images  are 
viewed  as  a  video  sequence,  these  elliptical  features  exhibit 
a  distinctive  outward  propagation  pattern  originating  from  a 
well-defined  center.  Unlike  prior  observations  of  such  con¬ 
centric  wave  structures  (Taylor  and  Hapgood,  1988;  Dewan 
et  al.,  1998),  the  concentric  wave  patterns  here  are  unusual 
in  that  they  form  complete  360°  closed  rings.  The  reason  for 
this  surprising  difference  will  become  clear  in  Section  5.2.3. 

The  gravity  wave  structure  shown  in  Fig.  3  may  be 
mapped  onto  a  fiat  surface  to  determine  the  relationship  of 
the  various  features  in  the  figure  to  the  underlying  thun¬ 
derstorm.  We  projected  each  pixel  in  the  image  onto  the 
85  km  tangent  plane  assuming  the  analysis  strip  (see  Sec¬ 
tion  5.2.1  below)  lies  in  the  tangent  plane.  Fig.  4  shows 
the  resultant  map  from  Fig.  3  laid  over  the  corresponding 
GOES  weather  map  of  the  Nebraska  storm.  This  projection 
clearly  associates  the  gravity  wave  modulation  structures 
with  the  underlying  thunderstorm. 


5.  Analysis 

5.1.  Identification  of  OH  Meinel  emissions 

The  video  sequence  shows  circular  waves  radiating 
outward  from  a  region  centered  over  the  thunderstorm  in 
eastern  Nebraska.  We  interpret  these  propagating  waves 
as  being  mesospheric  nightglow  emissions  modulated  by 
gravity  waves  generated  by  the  underlying  thunderstorm. 


There  are  several  nightglow  layers  in  the  mesosphere 
that  are  candidate  sources  for  the  observed  emissions. 
Only  two  fall  within  our  pass  band.  The  700-900  nm 
OH(X2F1,Av  =  3,4,5)  Meinel  sequences  indicated  in  Fig. 
1  emit  from  a  thin  layer  centered  at  an  altitude  of  ~85  km. 
The  O2  Atmospheric  Bands  02(hlXg  — ►  XlL~  )(0  -  0, 1 ) 
fall  within  our  pass  band,  but  the  atmosphere  is  opti¬ 
cally  thick  for  the  (0,0)  component  of  these  emissions  at 
761.9  nm  and  it  is  therefore  not  observed  at  the  ground. 
The  brightness  of  the  (0, 1 )  emissions  at  859.8  nm  is 
typically  ~300  R  (Slanger  et  al.,  2000),  and  therefore 
substantially  weaker  than  the  intensities  7~1  k R  detected 
here  (see  Section  5.4.1  below).  The  remaining  two  night¬ 
glow  sources,  the  Infrared  Atmospheric  Bands  02(al  Ag  — ► 
Ar3Z‘“,  At;  =  1)  -f  hv(  1.5  pm)  and  the  neutral  sodium  Na 
(2P?/2,3/2  ~ ^  1Si/2)  +  hv( 589  nm)  D-line  emissions,  fall 
outside  the  instrument  pass  band.  We  conclude  that  the  ob¬ 
served  emissions  are  most  likely  to  be  predominantly  from 
the  OH  Meinel  bands. 

5.2.  Characterization  of  gravity  wave  parameters 

5.2.1.  Analysis  method 

To  determine  the  phase  speed  and  wavelength  of  the  out¬ 
wardly  expanding  ripples,  we  established  a  horizontal  anal¬ 
ysis  strip  across  the  center  of  the  ripple  pattern  in  the  image, 
indicated  by  a  set  of  horizontal  dashed  lines  passing  through 
the  center  of  the  circular  ripple  patterns  in  Fig.  3.  For  each 
column  the  intensity  was  summed  across  4  rows  of  the  strip 
to  produce  a  horizontal  emission  profile.  By  assuming  the 
emissions  originate  at  an  altitude  of  85  km  corresponding  to 
the  OH  layer,  this  could  be  converted  to  the  corresponding 
horizontal  emission  profile  as  a  function  of  transverse  dis¬ 
tance  in  km.  One  such  profile  was  computed  for  each  30  s 
image  during  the  observing  interval. 

Fig.  5  shows  the  results  of  stacking  horizontal  emis¬ 
sion  profiles  from  successive  images  into  a  color-coded 
position-time  plot,  where  the  horizontal  scale  is  the  dis¬ 
tance  in  km  from  the  left  edge  of  the  image  at  the  position 
of  the  dashed  lines  in  Fig.  3.  Time  runs  vertically  from 
bottom  to  top  beginning  at  0400  UT.  The  outward  propa¬ 
gation  of  the  waves  from  the  central  region  of  the  figure 
with  increasing  time  is  clearly  discernible  in  this  plot. 
The  short  upward  streaks  in  the  plot  are  stars  as  they  pass 
through  the  (fixed)  analysis  strip.  Sprites  intersecting  the 
analysis  strip  produce  signatures  of  short  horizontal  strips 
of  width  equal  the  width  of  the  sprite.  However,  not  ev¬ 
ery  sprite  that  was  captured  in  the  images  intersected  the 
analysis  strip,  so  the  number  of  sprite  signatures  seen  here 
undercounts  the  actual  number  observed  of  approximately 
twenty. 

5.2.2.  Wave  period  and  horizontal  wavelength 

The  wave  period  and  horizontal  wavelength  may  be  read 
directly  from  Fig.  5,  and  are  approximately  t%8-1  1  min 
and  Ah%40  km,  respectively,  at  1  h  into  the  plot.  The 
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Fig.  5.  Color-coded  position-time  plot  of  emission  intensity  along  the  horizontal  analysis  strip  indicated  in  Fig.  2.  The  color  axis  runs 
through  purple-red-yellow  with  increasing  intensity.  The  waves  appear  to  propagate  outward  and  upward  from  the  central  vertical  region 
of  the  figure.  The  wavelengths  and  periods  of  the  waves  are  readily  discernible  from  this  plot,  yielding  a  phase  speed  of  about  60  m/s. 
The  right-to-left  drift  of  the  bright  (yellow)  central  region  with  increasing  time  tracks  the  movement  of  the  underlying  thunderstorm  across 
Nebraska.  Stars  and  sprites  that  appear  within  the  analysis  strip  are  labeled. 


corresponding  phase  speed  is  60-70  m/s.  This  period  and 
wavelength  differ  from  the  periods  of  5-13  min  and  wave¬ 
lengths  of  14-25  km  previously  reported  by  Taylor  and  Hill 
( 1991 ),  but  the  phase  speed  is  comparable.  In  Section  5.2.4 
below  we  analyze  the  time  evolution  of  the  propagation 
parameters  that  is  evident  in  Fig.  5. 

5.2.3.  Vertical  wavelength 

The  horizontal  component  cPh  of  phase  velocity  cp  is 
measured  to  be  cph^60  m/s  and  horizontal  wavelength 
4~40  km.  These  results  assume  emissions  from  an  altitude 
of  85  km.  The  resulting  period  t  is  670  s,  or  1 1  min. 

From  this  information  and  the  dispersion  relation  for  grav¬ 
ity  waves,  one  can  infer  other  characteristics  of  the  waves 
including  the  vertical  wavelength  kz.  Gravity  waves  orig¬ 
inating  in  the  lower  atmosphere  propagate  upward  with  a 
group  velocity  cg  which  is  orthogonal  to  the  phase  veloc¬ 
ity  cp,  or  the  wave  vector  k  =  (£h,fc)  =  (2n/k^2n/kz)y 
where  4  and  kz  are  the  horizontal  and  vertical  wave  num¬ 
bers,  and  4  and  4  are  the  corresponding  wavelengths. 
Hence,  cg  is  directed  along  the  phase  fronts,  and  the  angle 


of  elevation  (p  of  the  group  velocity  is  given  by 


tan  tp  =  — 


kz' 


(3) 


Short-period  waves  for  which  cph  <cSy  the  speed  of  sound, 
satisfy  the  dispersion  relation 


a j2 


,  42  2 

-  1  =  JS  =  COt  <P> 


(4) 


where  N  is  the  Brunt-Vaisala  frequency  and  co  =  2n/x  is  the 
gravity  wave  frequency.  It  follows  that: 

_2  \  “*l/2 


tan 


-OK 


(5) 


Assuming  a  Brunt  period  tb  =  2n/N&5  min  and  using  Eq. 
(3),  expression  (5)  yields 

tan  (p  =  ^  =  0.51,  (6) 

4 


with  the  result  that  tp  —  27°  and  4  =  20  km. 

Assuming  the  emissions  originate  in  the  OH  Meinel 
bands,  the  vertical  thickness  H  of  the  emission  layer  is 
roughly  8  km  (cf.  Makhlouf  et  al.,  1995).  Hence,  the  phase 
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of  the  wave  changes  very  little  in  the  z  direction  across  the 
layer  ( H  <  xr/2),  and  one  would  expect  to  have  little  “ob¬ 
servational  filtering”  of  the  wave  modulation  by  destructive 
interference  or  cancellation  effects  of  emissions  originat¬ 
ing  from  different  portions  of  the  emitting  layer.  This  is 
consistent  with  the  observed  gross  spatial  structure  of  the 
gravity  wave  emissions,  where  there  is  little  variation  in 
the  modulation  amplitude  across  the  entire  structure.  The 
large  vertical  wavelength  compared  to  H  is  the  reason  why 
the  concentric  wave  fronts  form  complete  closed  patterns, 
unlike  the  prior  observations.  By  way  of  contrast,  if  the 
vertical  wavelength  Az  is  less  than  //,  the  portions  of  the 
pattern  where  the  line-of-sight  was  aligned  more  closely 
with  the  wave  fronts  would  be  brighter  (Alexander  and 
Picard,  1999).  Such  a  thunderstorm-generated  wave  pattern 
with  varying  contrast  can  be  seen  in  the  MSX  satellite 
observations  of  stratospheric  gravity  waves  in  4.3-pm  CO2 
emission  (Dewan  et  al.,  1998).  There,  however,  the  effec¬ 
tive  thickness  of  the  emitting  layer  is  much  greater  (15- 
20  km)  (Picard  et  al.,  1998). 

5.2.4.  Time  evolution 

Over  the  observation  interval  the  horizontal  phase  speed 
of  the  waves  very  clearly  diminishes,  as  shown  in  Fig.  5.  At 
times  near  70  min  the  slopes  marked  in  this  figure  indicate 
a  phase  velocity  cp h  of  about  60  m/s,  corresponding  to  a 
horizontal  wavelength  Ah~40  km  and  a  period  r  =  1 1  min. 
However,  at  the  bottom  of  this  figure,  near  the  beginning  of 
the  analysis  interval  at  ~10  min  or  less,  the  slope  is  flatter, 
indicating  a  phase  velocity  of  about  85  m/s  and  a  horizon¬ 
tal  wavelength  of  about  50  km.  Using  r  =  Ah/cph  we  obtain 
t«9.8  min;  thus  both  wave  frequency  and  wavelength  were 
higher  in  the  earlier  part  of  the  storm.  Using  Eqs.  (5)  and 
(6)  we  can  calculate  Az  from  the  values  of  Ah  and  cPh.  For 
the  early  time  we  find  Az  =  30  km.  This  is  in  contrast  to  the 
value  Az  =  20  km  that  is  found  for  the  later  time.  In  Table  2 
we  summarize  these  parameters  and  compare  them  with  pre¬ 
vious  results  reported  by  Taylor  and  Hapgood  (1988)  and 
results  obtained  from  MSX  observations  on  two  different 
days  (Dewan  et  al.,  1998). 

Wavelength  parameters  have  been  related  to  the  physical 
properties  of  thunderstorms  by  Alexander  et  al.  (1995)  and 
Pandya  and  Alexander  (1999).  The  period  r  is  the  period  of 
mechanical  oscillation  of  the  thunderstorm,  i.e.  the  period 
of  the  oscillating  convective  updrafts  and  downdrafts  within 
the  storm  that  impinge  upon  the  tropopause.  The  value  of  Az 
was  shown  by  Salby  and  Garcia  (1987)  to  be  determined  by 
the  characteristic  vertical  dimension  of  ~  1 0  km  thickness  of 
the  thunderstorm’s  thermal  forcing  layer  in  the  troposphere. 
In  fact  the  vertical  extent  of  the  storm  should  correspond  to 
a  half  wavelength  A-/2,  but  in  propagating  from  the  tropo¬ 
sphere  to  the  stratosphere  the  vertical  wavelength  shortens 
by  about  a  factor  of  2  due  to  wave  refraction  (Alexander 
et  al.,  1995),  so  that  Az  in  the  stratosphere  and  above  is 
on  the  order  of  the  vertical  extent  of  the  storm.  On  the  basis 
of  these  considerations  we  see  from  Table  2  that  the 


oscillation  frequency  and  vertical  extent  of  the  storm 
both  decrease  in  time  during  the  life  of  the  storm.  The 
values  of  20  and  30  km  inferred  from  the  observations 
appear  to  be  too  large  to  represent  actual  storm  extents, 
but  the  expected  qualitative  evolution  of  the  thunderstorm 
parameters,  consisting  of  slowing  of  the  mechanical  os¬ 
cillation  and  shrinking  of  the  storm’s  vertical  dimension, 
seems  to  be  firmly  established  by  the  present  measure¬ 
ments.  Additional  wavelength  stretching  due  to  wave 
refraction  occurring  when  the  wave  propagates  into  the 
less  stable  mesosphere  may  account  for  a  portion  of  the 
discrepancy. 

5.3.  Other  wave  characteristics 

From  the  wave  propagation  elevation  angle  <p ,  the  ground 
range  between  the  source  and  the  region  where  the  wave 
pattern  is  observed  in  the  nightglow  may  be  determined.  If 
we  assume  rectilinear  propagation  without  refraction  from 
wind  shears,  for  a  thunderstorm  source  altitude  of  15  km 
near  the  tropopause  and  an  OH  emission  altitude  of  85  km 
the  ground  range  is  70  km/0.5 1«  140  km.  As  can  be  seen  in 
Fig.  3  this  corresponds  roughly  to  the  radius  of  the  observed 
wave  pattern,  so  the  characteristic  horizontal  scale  size  of 
the  gravity  wave  effects  is  on  the  order  of  a  few  times  the 
height  of  the  nightglow  layer  above  the  thunderstorm  source. 

The  transit  time  from  the  source  to  the  altitude  where 
the  effects  are  observed  may  be  estimated  to  be  approxi¬ 
mately  70  km/cgz,  where  cgz  is  the  vertical  component  of 
the  group  velocity  of  the  gravity  wave,  given  by  cg2  = 
oco/dkz  =  (Az/t)cos2  </>,  where  co  =  2n/x  is  the  wave  fre¬ 
quency  (Dewan  et  al.,  1998).  Substituting  the  current  values 
of  the  wave  parameters,  we  calculate  the  vertical  transit  time 
from  the  thunderstorm  to  the  mesopause  OH  layer  to  be  ap¬ 
proximately  40-50  min.  Hence,  in  comparing  thunderstorm 
characteristics  and  gravity  wave  parameters,  storms  that  can 
act  as  sources  should  be  determined  by  looking  at  GOES 
images  40-50  min  before  the  wave  observation. 

Over  the  observation  interval  the  underlying  thunderstorm 
moved  at  a  speed  of  approximately  15-20  m/s,  roughly 
0.2- 0.3  the  horizontal  phase  speed  of  the  gravity  waves. 
The  region  of  the  nightglow  layer  lying  directly  above  the 
thunderstorm  source  moves  with  the  thunderstorm,  which 
from  Fig.  3  is  seen  to  be  predominantly  from  right-to-left  as 
observed  from  Bear  Mt.  looking  southeastward.  This  move¬ 
ment  of  the  thunderstorm  gravity  source  is  reflected  in  the 
nightglow  modulations,  which  can  be  seen  in  Fig.  5  as  a 
right-to-left  drift  in  the  brightest  portion  of  the  emissions  as 
time  progresses  from  bottom-to-top. 

5.4.  OH  brightness  perturbations  and  associated 
temperature  perturbations 

5.4.1.  Gravity  waves 

Absolute  intensity  calibrations  were  not  available  for 
the  Photometries  camera,  but  the  intensities  of  the  OH 
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Table  2 


Comparison  of  reported  gravity  wave  parameters 


4  (km) 

A-  (km) 

t  (min) 

cph  (m/s) 

Cgh  (m/s) 

cgz  (m/s) 

At  (min) 

<P  (deg) 

Before  0430 

50 

30 

10 

85 

58 

27 

43 

31 

After  0500 

40 

20 

11 

61 

48 

24 

49 

27 

MSX  11/13 

25 

16 

9 

46 

32 

21 

55 

33 

MSX  10/2 

50 

18 

15 

56 

50 

18 

64 

20 

T&  H  1988 

26 

7 

17 

38 

23 

11 

105 

15 

Ah:  Horizontal  wavelength. 

Ar:  Vertical  wavelength, 
i:  Wave  period. 

cph;  Wave  horizontal  phase  speed. 

Cgh :  Wave  horizontal  group  speed. 
cgz :  Wave  vertical  group  speed. 

At:  Wave-energy  propagation  time  from  15  to  85  km  altitude. 
q>:  Elevation  angle  of  group-velocity  vector. 


emissions  may  be  estimated  from  the  brightness  of  the 
sprites  simultaneously  recorded  with  the  broad  band,  high 
speed  imager  operated  during  the  Sprites99  campaign 
(Stenbaek-Nielsen  et  al.,  2000).  Sprites  typically  persist  for 
a  few  (~3-10  ms)  and  strongly  saturate  the  3  MR  limit  of 
the  ICCD.  The  large  sprites  observed  during  this  storm  are 
estimated  to  possess  brightness  ~10  MR.  The  images  ob¬ 
tained  with  the  Photometries  NIR  CCD  register  the  apparent 
nightglow  brightness  to  be  approximately  one-half  that  of 
the  sprites.  By  taking  into  account  the  sprite  duration  (~5  ms 
is  a  typical  value)  and  the  image  integration  time  (25  s), 
the  brightness  of  the  gravity  wave  modulated  OH  emis¬ 
sions  is  estimated  to  be  (±)(5  ms/25  s)(  10  MR)~1  kR. 
This  brightness  is  comparable  to  the  OH  brightness  of 
1-2  kR  over  a  similar  wavelength  pass  band  reported 
by  Yee  et  al.  (1991)  for  observations  made  during  the 
ALOHA-90  campaign.  We  shall  henceforth  use  1  kR 
for  the  mean  brightness  of  the  OH  nightglow  in  our 
observations,  recognizing  that  this  is  an  approximate 
figure. 

The  optical  modulation  amplitude  of  the  outwardly  ex¬ 
panding  waves  was  measured  to  be  approximately  30%  peak 
to  peak  compared  to  the  background  mean  emission  inten¬ 
sity.  From  Eq.  (2)  and  using  values  of  the  phase  speed 
given  in  Table  2,  we  estimate  rj^l.  The  resultant  tempera¬ 
ture  perturbation  associated  with  the  concentrically  expand¬ 
ing  waves  is  ~4%  of  the  ambient,  or  ~7  K. 


5.4.2.  Inferred  upper  limit  for  energy  deposited  by  sprites 
in  the  OH  layer 

The  expected  emission  signature  of  a  localized  transient 
temperature  perturbation  in  the  nightglow  layer  that  sig¬ 
nificantly  alters  the  OH  emissivity  is  a  pulse  that  spreads 
radially  outward  from  the  source  at  the  acoustic  speed.  We 
examined  the  observations  carefully  for  the  occurrence  of 


small,  radially  expanding  ripples  originating  from  the  point 
where  sprites  penetrated  the  OH  nightglow  layer.  In  Fig. 
5  they  would  be  expected  to  appear  as  upward-slanting 
“wings”  originating  from  each  sprite  that  fell  within  the 
analysis  strip.  No  sprite  signatures  of  this  type  were  de¬ 
tected  in  the  OH  emissions  at  the  level  of  2%  of  the  back¬ 
ground  at  that  location.  Using  an  assumed  OH  brightness  of 
~  1  kR  this  negative  result  permits  us  to  establish  an  obser¬ 
vational  upper  limit  of  ~20  R  for  the  brightness  perturbation 
in  the  OH  produced  by  a  large  sprite  at  altitudes  ~85  km. 
From  Eq.  (2)  et  seq.  and  using  values  of  the  wave  period 
and  phase  speed  given  in  Table  2,  we  estimate  r\^l.  The 
corresponding  brightness  perturbation  AI/I  <;  2%  trans¬ 
lates  to  an  upper  limit  for  the  corresponding  neutral  temper¬ 
ature  perturbation  AT/T  <  0.003,  or  A  T  <  0.5  K.  at  the 
altitude  of  the  mesopause. 

The  upper  limit  AT <0.5  K  for  the  temperature 
perturbation  may  be  converted  to  an  upper  limit  for  the 
thermal  energy  deposited  into  the  local  medium  by  the 
sprite.  If  we  assume  that  most  of  the  energy  associated 
with  electron  impact  excitation  is  ultimately  transferred 
into  thermal  energy  by  collisions  with  neutrals,  then  the 
total  energy  associated  with  this  process  may  be  estimated 
as  AE^NnksVAT,  where  Nn  is  the  number  density  of 
neutrals,  k&  =  1.38  x  10-23  J/K  is  Boltzmann’s  constant, 
and  V  is  the  effective  volume  of  the  portion  of  the  sprite 
intersecting  the  OH  layer.  This  volume  is  estimated  using 
the  vertical  thickness  ~7  km  of  the  OH  emission  layer  and 
a  horizontal  cross  section  of  108  m2  of  the  sprite  where  it 
penetrates  the  OH  layer,  yielding  K~7  x  1011  m3.  With 
neutral  density  Nn  =  2  x  1020/m3  at  an  altitude  of  85  km 
and  AT <0.5  K,  we  arrive  at  an  upper  limit  A£<1  GJ 
for  the  associated  total  thermal  energy  deposited  in  the 
OH  layer.  As  will  be  shown  in  the  next  section,  this  upper 
limit  is  very  much  greater  than  the  characteristic  energies 
deposited  by  sprites  in  this  region. 
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5.4.3.  Estimate  of  total  energy  deposited  by  sprites  into 
the  mesosphere 

To  estimate  the  total  energy  deposited  by  the  large 
sprites  observed  in  association  with  the  thunderstorm 
studied  here  we  use  the  method  presented  by  Heavner  et 
al.  (2000),  with  parameters  adapted  to  the  present  prob¬ 
lem.  The  method  derives  from  estimates  of  total  optical 
energy  and  electron  temperature,  and  scales  the  problem 
to  include  contributions  from  unobserved  infrared  emis¬ 
sions,  quenching,  and  inelastic  collision  effects  from  N2. 
The  optical  energy  in  the  present  case  may  be  estimated 
from  the  total  brightness  and  cross-section  of  the  sprite, 
and  the  duration  of  the  optical  emissions.  The  bright¬ 
ness  B( R)  in  Rayleighs  is  given  by  (Chamberlain,  1961, 
Appendix  II)  Z?(R)  =  10-6  /path  F(r)dr,  where  the  integral 
is  taken  along  the  line  of  sight  through  the  source,  and 
where  F(r)  is  the  volume  emission  rate.  For  a  brightness  of 
10  MR,  cross-section  of  (10  km)2,  and  duration  3  ms  we 
obtain  a  total  number  of  photons  3  x  1022  emitted  by  the 
sprite  in  the  imager  pass  band.  Assuming  an  average  photon 
wavelength  —800  nm— 1.5  eV  for  the  observed  photons, 
the  total  optical  energy  emitted  is  4.5  x  1022  eV— 7.5  kJ. 
Using  the  most  recent  estimates  for  the  average  electron 
energy  —2  eV  (Morrill  et  al.,  2002)  and  N2  vibrational  dis¬ 
tributions  (Bucsela  et  al.,  2003),  this  results  in  an  estimated 
upper  limit  of  —3.3  MJ  for  the  total  energy.  It  should  be 
noted  that  this  method  of  determining  total  energy  based  on 
optical  brightness  is  very  sensitively  dependent  on  electron 
temperature  when  it  is  much  less  than  2  eV  and  should  be 
used  with  caution  in  this  regime. 

The  estimate  of  —3.3  MJ  for  the  total  energy  compares 
favorably  with  the  total  energy  of  —2  MJ  independently 
estimated  by  Armstrong  et  al.  (2001 )  based  on  a  comprehen¬ 
sive  model  calculation  for  sprites  of  similar  size  and  bright¬ 
ness.  It  is  also  of  the  same  order  of  magnitude  as  the  total 
electrostatic  energy  available  in  the  transient  electric  field 
following  a  large  lightning  discharge.  This  latter  may  be  es¬ 
timated  from  W  =  sqE2  F/2,  where  eo  =  8.86  x  10-12  F/m 
is  the  permittivity  of  free  space,  E  is  the  amplitude  of  the 
transient  electric  field,  and  V  is  the  volume  of  the  sprite. 
Using  £—300  V/m  at  75  km  (Pasko  et  al.,  1997a)  and  as¬ 
suming  a  volume  V  =  1012  m3  we  obtain  IT- 800  kJ.  The 
range  of  values  800  kJ-3.3  MJ  reflects  differences  in  model 
assumptions,  as  well  as  uncertainties  in  the  brightness/size 
distributions  of  sprites,  but  are  consistent  enough  at  the  level 
of  order  of  magnitude  to  suggest  a  characteristic  energy  of 
~1  MJ  for  large,  bright  sprites.  The  uncertainty  factor  as¬ 
sociated  with  this  value  remains  substantial  owing  to  the 
wide  variety  of  assumptions  used  in  the  estimates,  so  the  es¬ 
timate  is  probably  no  better  than  order  of  magnitude.  How¬ 
ever,  even  taking  into  account  model  uncertainties  this  is 
significantly  below  previous  upper  limit  estimates  of  —  1  GJ 
total  energy  in  sprites  (e.g.,  Heavner  et  al.,  2000;  Sentman 
et  al.,  2002),  so  characteristic  energy  estimates  should  be  re¬ 
vised  downward  to  the  current  order  of  magnitude  estimate 
of -1-10  MJ. 


6.  Discussion 

The  characteristic  total  energy  — 1-10  MJ  of  a  sprite  es¬ 
timated  above  is  clearly  much  smaller  than  the  estimated 
upper  limit  —  1  GJ  needed  to  produce  a  detectable  response 
in  the  OH  brightness  using  the  reactions  ( 1 )  and  (I').  If  we 
further  assume  that  local  dissipation  of  sprite  energy  takes 
place  near  the  —75  km  altitude  region  of  brightest  optical 
emissions,  then  only  some  fraction  of  the  total  energy  will 
be  deposited  in  the  higher  OH  layer  at  —85  km,  thus  fur¬ 
ther  widening  this  difference.  The  size  of  the  discrepancy 
between  the  required  and  available  energy  appears  to  effec¬ 
tively  rule  out  electrical  effects  associated  with  sprites  play¬ 
ing  a  significant  role  in  OH  airglow  processes  by  means  of 
localized  thermal  heating  alone. 

One  may  speculate  that  other  mechanisms  could  produce 
a  signature.  Here,  we  have  assumed  the  source  mechanism 
for  impulsive  sprite  production  of  acoustic  waves  is  local¬ 
ized  Joule  heating  within  the  sprite  volume.  The  expected 
effect  was  the  launching  of  transient  pulses  or  wave  pack¬ 
ets  from  the  sprite  source  where  it  penetrated  the  OH  air- 
glow  layer,  and  the  expected  signature  was  modulation  of 
the  temperature  and  the  ambient  OH  emissions  by  the  wave 
packets.  Joule  heating  of  the  neutrals  by  quasi-electrostatic 
fields  or  induced  currents  in  the  sprite  volume  is  the  most 
likely  mechanism  for  inducing  localized  bulk  changes  to  the 
medium.  However,  other  consequential  effects  might  include 
perturbation  to  the  airglow  chemistry  of  nocturnal  bands  be¬ 
sides  OH  that  lie  within  the  mesopause  vicinity  of  sprites, 
such  as  sodium  or  the  oxygen  bands,  although  these  would 
be  difficult  to  study  from  the  ground. 

Other  signatures  might  include  local  spatial/temporal  ef¬ 
fects  from  the  medium  itself.  We  have  focused  here  on 
the  gravity  wave  signature,  but  one  might  also  consider  the 
infrasound  branch  of  the  acoustic  dispersion  relation  and 
search  for  modulations  in  the  optical  emissions  at  infra- 
sound  frequencies  (—0.1-1  Hz)  using  ground  based  micro¬ 
barograph  measurements  or  HF  measurements  of  the  lower 
ionosphere. 

7.  Conclusions 

7.1.  Principal  results 

We  have  performed  a  simplified  analysis  on  the  OH  grav¬ 
ity  wave  structures  observed  simultaneously  with  sprites  to 
determine  the  basic  wavelength  and  period  parameters  of  the 
gravity  waves,  and  to  extract  the  corresponding  temperature 
perturbations  implied  by  these  parameters.  The  absence  of 
signatures  in  OH  that  could  be  associated  with  sprites  was 
used  to  estimate  an  upper  limit  on  the  energy  perturbation 
associated  with  the  sprite.  Separate  calculation  of  the  energy 
in  a  sprite  based  on  optical  emission  intensities  showed  that 
the  energy  deposited  by  a  sprite  in  the  mesosphere  is  well 
below  that  required  to  produce  a  detectable  perturbation  in 
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the  OH  emissions.  Our  principle  results  are: 

1.  Sprites  and  circular,  outwardly  expanding  patterns  of 
optical  emissions  were  simultaneously  observed  in  NIR 
720-920  nm  over  an  intense  Nebraska  thunderstorm  on 
the  night  of  18  August  1999.  The  expanding  circular 
emission  patterns  are  interpreted  as  modulation  of  the 
nightglow  OH  (X2  77,  At;  =  3,4,5)  Meinel  sequences  at 
altitudes  of  ~85  km  by  gravity  waves  propagating  up¬ 
ward  from  the  underlying  thunderstorm. 

2.  The  horizontal  wavelength  and  period  of  the  gravity 
waves  underwent  an  evolution  during  the  2  h  observ¬ 
ing  interval.  During  the  early  stages  of  the  storm  the 
horizontal  wavelength  was  measured  to  be  50  km  and 
the  period  was  9.8  min,  corresponding  to  a  phase  speed 
of  85  m/s.  The  vertical  wavelength  was  determined  to 
be  30  km.  Approximately  1  h  later  these  parameters 
had  evolved  to  40  km  wavelength,  1 1  min  period,  a 
markedly  slower  phase  speed  of  60  m/s  and  a  smaller 
vertical  wavelength  of  20  km.  The  large  vertical  wave¬ 
length  compared  to  the  7-km  thickness  of  the  OH 
emission  layer  results  in  an  unusual  complete  concentric 
wave  pattern  being  observed,  free  of  observational  filter¬ 
ing  effects.  The  vertical  wavelength  of  the  gravity  waves 
is  several  times  the  estimated  vertical  dimension  of  the 
convective  pump  of  the  underlying  thunderstorm  source, 
but  the  qualitative  temporal  change  of  wave  parameters 
is  in  accord  with  present  understanding  of  thunderstorm 
dynamics. 

3.  The  average  brightness  of  the  OH  layer  was  estimated 
from  comparison  with  simultaneously  observed  sprites 
to  be  approximately  1  kR.  The  peak-to-peak  modulation 
intensity  of  the  gravity  wave  modulated  OH  emissions 
was  approximately  30%  of  the  mean  during  the  most  in¬ 
tense  portion  of  the  event.  This  level  of  modulation  cor¬ 
responds  to  an  underlying  neutral  pressure/temperature 
perturbation  of  approximately  4%  about  the  mean.  For  a 
mesopause  temperature  of  ~170  K  the  associated  tem¬ 
perature  perturbation  is  ~7  K. 

4.  Sprites  observed  coincident  with  the  gravity  wave  effects 
produced  no  discernible  perturbations  or  distinctive  sig¬ 
natures  in  the  OH  emissions  at  the  2%  brightness  level. 
This  corresponds  to  an  upper  limit  in  the  brightness  per¬ 
turbation  of  ~20  R,  and  an  upper  limit  for  sprite  heating 
of  AT  <  0.5  K  at  the  altitude  of  the  mesopause.  The  cor¬ 
responding  total  thermal  energy  deposition  is  A E  <  1  GJ 
if  uniform  emission  brightness  is  assumed  within  the 
volume  of  a  sprite.  This  upper  limit  on  the  energy  lies 
well  above  the  total  estimated  energy  of  ~  1—10  MJ 
deposited  in  the  mesosphere  by  a  sprite.  The  wide 
difference  between  the  ~1— 10  MJ  characteristic  total 
energy  of  bright  sprites  and  the  ~  1  G J  energy  required 
to  produce  a  detectable  optical  response  in  OH  emis¬ 
sions  indicates  that  sprites  do  not  produce  a  significant 
effect  on  mesospheric  OH  processes  by  means  of  thermal 
heating. 


7.2.  Further  comments  on  gravity  wave  observations 

While  the  principal  motivation  for  undertaking  the  present 
work  was  to  use  possible  gravity  wave  signatures  in  air- 
glow  as  a  tool  to  extract  the  energy  deposition  of  sprites 
in  the  mesosphere,  the  image  data  set  provides  a  valuable 
resource  in  its  own  right  for  studying  the  gravity  waves 
associated  with  the  storm  of  18  August  1999.  Additional 
analyses  that  are  possible  with  this  data  include  estimat¬ 
ing  the  power  in  the  waves  and  comparing  it  with  what 
is  available  in  the  storm,  considering  the  efficiency  of  the 
coupling  to  the  atmosphere.  This  requires  taking  into  ac¬ 
count  the  wave-generation  mechanisms.  The  three  princi¬ 
pal  mechanisms  by  which  storms  can  generate  waves  are 
[M.J.  Alexander,  private  communication]: 

(a)  The  “ moving  mountain ”:  Shear  flow  over  a  dome  of 
air  associated  with  a  moving  convective  cell,  that  is  a  “mov¬ 
ing  mountain”  figuratively,  can  launch  “mountain”  waves 
upward. 

(b)  Transient  pumping :  This  is  the  mechanical  forcing 
effect  on  the  atmosphere  from  the  body  force  associated  with 
the  moving  expanding  air  in  the  cell.  There  are  two  limits. 
If  the  forcing  is  impulsive,  then  it  contains  all  frequen¬ 
cies  and  will  tend  to  generate  a  broad  spectrum  of  waves. 
According  to  the  dispersion  relation,  the  waves  will  be  seg¬ 
regated  in  direction  by  their  frequency,  with  short-period 
waves  propagating  more  nearly  vertically  and  longer  pe¬ 
riod  waves  propagating  in  a  more  horizontal  direction. 
On  the  other  hand,  there  is  evidence  that  quasi-sinusoidal 
oscillations  of  the  convective  air  column  occur  in  thunder¬ 
storms.  Such  oscillations  would  give  rise  to  periodic  forcing 
and  generate  quasi-sinusoidal  waves.  These  waves  would  be 
launched  at  an  angle  determined  by  the  period  of  the  forcing 
and  the  wave  dispersion  relation.  We  have  tacitly  assumed 
this  situation  above  in  Section  5.2,  as  well  as  in  Dewan 
et  al.  (1998). 

(c)  Transient  heat  source :  This  is  the  driving  term  in  the 
thermal  equation  from  the  heat  release  in  the  storm.  It  tends 
to  generate  waves  of  given  vertical  wavelength  2_-  rather  than 
given  frequency  or  period.  The  heated  air  column  generates 
waves  whose  half-wavelength  kz!2  is  equal  to  the  overall 
length  of  the  column,  as  discussed  in  Section  5.2.5. 

These  mechanisms  all  operate  simultaneously  with  vary¬ 
ing  relative  contributions  to  the  total  energy  input  into  the 
wave  field  and  are  all  seen  in  the  numerical  simulations 
of  convective  wave  generation  (Alexander  et  al.,  1995; 
Pandya  and  Alexander,  1999).  For  the  present  set  of  ob¬ 
servations  of  concentrically  expanding  waves  emanating 
from  a  point  above  a  thunderstorm,  the  original  explanation 
offered  by  Taylor  and  Hapgood  (1988)  appears  to  hold.  In 
that  report,  the  occurrence  of  concentric  circular  waves  in 
the  OH  airglow  assumed  that  the  thunderstorm  was  a  quasi- 
monochromatic  source  of  gravity  waves,  and  evidence  for 
this  was  cited  in  the  work  of  Anderson  (1960).  Most  re¬ 
cently,  Piani  et  al.  (2001 )  performed  three-dimensional  simu¬ 
lations  of  gravity  waves  generated  by  a  large  connective 
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thunderstorm  and  obtained  results  in  striking  agreement 
with  observations  of  Taylor  and  Hapgood  (1988),  Dewan 
et  al.  (1998),  and  the  present  work.  Thus,  the  cause  of  the 
observed  concentric  patterns  appears  to  be  well  understood. 

Finally,  the  duration  of  excitation  at  the  center  of  the  out¬ 
wardly  expanding  gravity  wave  ripples  may  be  interpreted 
as  being  caused  by  periodic  or  quasi-periodic  ringing  at 
the  tropopause  due  to  pumping  by  the  buoyant  air  column 
in  the  convective  cell  below.  The  resultant  pressure  pulse 
is  transmitted  through  the  stratosphere  and  mesosphere  at 
the  pump  frequency.  The  combination  of  spatially  concen¬ 
tric  localization  of  the  gravity  wave  signatures  and  their 
quasi-monochromatic  temporal  oscillations  would  seem  to 
suggest  a  compact  source  at  the  tropopause  with  a  charac¬ 
teristic  pump  period  reflecting  the  convective  column  dy¬ 
namics.  The  long  duration  of  the  gravity  wave  signatures 
suggests  that  the  source  must  have  been  at  least  narrow 
band/quasi-pcriodic. 
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